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There is a claim that a static charged black hole (Reissner-Nordstrom black hole) can be overcharged by 
absorbing a charged test particle. If it is true, it might give a counter example to the weak cosmic censorship 
conjecture, which states that spacetime singularities are never observed by a distant observer. However, so far 
the proposed process has only been analyzed within a test particle approximation. Here we claim that the back 
reaction effects of a charged particle cannot be neglected when judging whether the suggested process is really 
a counter example to the cosmic censorship conjecture or not. Furthermore, we argue that all the back reaction 
effects can be properly taken into account when we consider the trajectory of a particle on the border between 
the plunge and bounce orbits. In such marginal cases we find that the Reissner-Nordstrom black hole can never 
be overcharged via the absorption of a charged particle. Since all the plunge orbits are expected to have a higher 
energy than the marginal orbit, we conclude that there is no supporting evidence that indicates the violation of 
£ — ■ the cosmic censorship in the proposed overcharging process. 
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OJQ. I. INTRODUCTION 

' General relativity is the most successful theory of gravity and it has brought us deep a understanding of spacetime. Never- 
\ theless, when we evolve the Einstein equations with a well-posed initial condition, singularities, at which general relativity and 
all established theories lose their predictability, are known to form. However, in most cases singularities are hidden by event 
horizons as in the case of black holes and cannot be seen by a distant observer. It requires extreme fine-tuning of the initial state 
\^ , or unphysical equation of state to produce naked singularities. This statement is known as the cosmic censorship conjecture 
proposed by Penrose 01 . Despite the tremendous efforts to prove whether the cosmic censorship conjecture is a generic property 
of classical general relativity, it still remains an open question 10]. 

According to the uniqueness theorem H, all stationary asymptotically flat black holes in Einstein-Maxwell system are de- 
scribed by Kerr-Newman solutions, which are specified uniquely by the mass M, the charge Q and the angular momentum J 
satisfying 
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^ ; M 2 >Q 2 + (J/M) 2 . (1) 

When the equality is saturated, the black holes are called extremal, and a naked singularity appears when M 2 < Q 2 + (J/M) 2 . 
Naively, it seems possible to form a naked singularity by throwing matter into a black hole to increase its charge and angular 
momentum. If it were really possible, we would say that the black hole is overcharged or overspinned via matter absorption, and 
it would give a counterexample to the cosmic censorship conjecture. 

A number of previous works related to overcharging and overspinning a black hole support the cosmic censorship conjecture. 
A pioneering work by Wald proved that neither overcharging nor overspinning is possible when a test (charged) particle plunges 
into an extremal black hole, and hence a naked singularity cannot be produced |4)j. If the particle is carrying charge or angular 
momentum sufficient to overcharge (or overspin) the black hole, it is not captured by the black hole because of the electro- 
magnetic or centrifugal repulsion force. Motivated by Wald's analysis, there have been many analyses attempting to supersaturate 
an extremal black hole by capturing a particle and a wave packet of a classical field. JH]. All these analyses indicate that such 
processes at most sustain the extremal condition if the particle is allowed to fall into the black hole. 1 

Recently, however, an alternative viewpoint was raised, suggesting the possibilities of overcharging a Reissner-Nordstrom 
(RN) black hole (01, of overspinning a Kerr black hole ! , and of overspinning or overcharging a Kerr-Newmann black hole H . 
The point is to consider the process of particle absorption with the initial black hole being prepared slightly below the extremal 
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1 There is an exception for this statement. De Felice and Yu have pointed out that even an extreme RN black hole can turn into a naked singularity if an 
electrically neutral spinning particle is sent to an extreme RN black hole with appropriate initial conditions ((J. 
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limit. At a first glance, these works are indicating that the violation of cosmic censorship conjecture is possible. However, as 
already mentioned in Refs. 0] and [H], the back reaction effects due to the presence of a particle are not taken into account in 
these analyses. In fact, it was emphasized in Ref. [ 1 00 that the finite size effect of a charged particle and the contribution of the 
interaction energy between the black hole and the particle can be important enough to protect the cosmic censorship conjecture. 
In Ref. Hill , the loss of energy due to gravitational radiation in the spin-up process suggested in Ref. (H was evaluated. Although 
the overspinning of a Kerr black hole can be avoided for some parameter choices by taking into account the effect of this radiative 
loss alone, the possibility of overspinning still remains for the other parameter choices. 

Our main purpose of this paper is to clarify the role of back reaction effects in the gedanken experiment discussed in Ref. 0]. 
Namely, we discuss whether a RN black hole can be overcharged or not by the capturing of a charged particle, taking into 
account all possible back reaction effects. The basic idea is to concentrate on the case in which the particle is in the orbit at the 
border between plunge and recoil orbits, which we call the marginal orbit. When the particle takes the marginal orbit, it will 
experience an unstable equilibrium configuration at the separatrix. Thanks to this property of the marginal orbit, one can use an 
exact solution of the Einstein-Maxwell equations to read the back reaction effects with sufficient accuracy for this equilibrium 
configuration. We will show that the total energy of the system is always greater than the total charge. Furthermore, using the 
black hole perturbation technique, one can evaluate the radiative loss of energy as the particle falls into the black hole from this 
unstable equilibrium configuration. We will find that the energy loss through this process is always negligible small. As a result, 
we will conclude that a particle in the marginal orbit cannot overcharge the RN black hole. Naively, the orbits that can plunge 
into the black hole will have higher energy than the marginal orbit. Hence, our result will indicate that the cosmic censorship 
conjecture is protected even if we consider the process suggested in Ref. 0. 

The outline of this paper is as follows. In Sec. [II] we briefly review the overcharging process suggested in Ref. 0]. After 
explaining the marginal orbit in Sec. [TTTJ we review the basic properties of an exact solution that corresponds to an unstable 
stationary equilibrium configuration composed of two charged sources, and we will prove that the total energy is always greater 
than the total charge for this configuration. In Sec. [IV] we will show that the radiative energy loss is always suppressed as long 
as a particle falls from an equilibrium position to the black hole. We summarize the results in Sec.[V] adding discussion about 
the implication of our results. 

In this paper, we use the units in which G = c = 1, but we explicitly write k 2 = 8ttG in Sec.[lV]and AppendixlDl respecting 
the original notations in Ref. Il2tl . The sign convention of the metric as (—,+,+,+). We adopt the definition of Riemann tensor 
and Ricci tensor given by R p ajlv := d^Y p a — d v Y p a + Y p yY^ a — Y p x Y* a and i? M „ := R a mv - Fourier components of ^>(t) are 

defined by {2tt)~^ 2 e iut ^(t)dt. 



II. OVERCHARGING A REISSNER-NORDSTROM BLACK HOLE WITHOUT BACK REACTION EFFECTS 

We briefly review the basic idea of overcharging a RN black hole by the infall of a charged particle, giving the precise meaning 
of "overcharging", based on the discussion in Ref. 0]. In this section, we first focus on the test particle case in which we neglect 
the back reaction effects. In this approximation we have a broad range of configurations whose final state exceeds the extremal 
bound. Later, we discuss how the results are modified by taking into account the back reaction effects. 



A. Test particle case 



We consider a point particle with mass fj, and charge q radially falling toward a nearly extremal RN black hole with mass M 
and charge Q. By assumption, these parameters satisfy fi < q -C Q < M. The background metric and the vector potential of a 
RN black hole are given by 

ds 2 = -f(r)di 2 + -^-dr 2 + r 2 (d6 2 + sin 2 ^ 2 ) , 
f{r) 

A = -~, A r = A e =A 4> = 0, (2) 



where 



/W:=1 _^ + «i (3) 



r 



Due to the spherical symmetry of a RN black hole, we can choose the coordinate system such that the trajectory of a particle 
is along the axis without loss of generality. Then, the trajectory is represented by z a (s) — (T(s), R(s),0, 0) with the proper 
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time s along the world line of the particle. The equations of motion of the particle are given by 

dR\ 2 1 /_ qQ x 



ds ) nf(R) V R 

The energy of the particle E is defined by E := —(dt) a (fiu a + qA a ) where (d t ) is the Killing field associated with the time 
coordinate t and u a := dz a /ds is the four velocity of the particle. Using the f-component of the equations of motion (0, it is 
shown that E is constant along the orbit. 

Radial orbits toward the black hole can be classified into three classes: plunge, bounce and marginal orbits. In the case of 
plunge orbits, a particle falls into the black hole horizon without changing the direction of motion. In the case of bounce orbits, 
there exists a turning point outside the horizon, at which (dR/ds) 2 = 0. The particle is reflected back at this point due to 
the electric repulsion force. The marginal orbit is on a separatrix between plunge and bounce orbits. In this case, the particle 
gradually approaches an unstable equilibrium position at r = ro- The position of ro and the energy of the marginal orbit are 
simultaneously determined by the conditions dR/ds = d 2 R/ds 2 = 0. 

In the above setup, if the following two conditions are satisfied, we would say that a RN black hole is possibly overcharged. 
The first condition is that the particle is in a plunge orbit. For plunge orbits there is no turning point where the radial velocity 
becomes zero. Therefore the condition is 

dR \ 2 

— ) = V (o) (R) > 0, for R> r + , (absorption condition) (6) 
ds J 



where r + := M + y M 2 — Q 2 is the radial coordinate of the event horizon. The second condition is on the total energy of the 
final state. After the absorption of the particle, the system will approach another RN geometry with mass Mfi na i = M + E and 
charge Qfi na i = Q + q. The condition that the final RN geometry exceeds the extremal bound is given by 

M + E < Q + q . (overcharging condition) (7) 

In Ref. 10], it was demonstrated that radial orbits in a rather wide range of parameter space satisfy the conditions (0 and (0. 
In fact, both the absorption and overcharging conditions are satisfied for 



1< b < a , c< \/a 2 - b 2 , (8) 

under the parametrization 

M := 1 + 2e 2 , Q := 1 , E := ae - 2be 2 , q := ae , fi := ce , (9) 

where < e< 1 and the coefficients a, b and c are assumed to be 0(1) real numbers. Here the overcharging is discussed at 
the level of 0(e 2 ), while the mass and the charge of a point particle are 0(eV Although the back reaction effects are completely 
neglected, the analysis in Ref. [Q] is in sharp contrast with previous works J5|. The previous analyses are restricted to O(e) and 
only the extremal limit is considered for the initial black hole. 



B. Back reaction effects on the overcharging process 

As was already emphasized in Ref. 10], the analysis in the previous section is not sufficient to conclude that a nearly extremal 
RN black hole can be overcharged through the absorption of a charged particle. We have to compute the total amount of energy 
in the final state of the whole system to the accuracy of 0(e 2 ), taking into account the effects of back reaction ifTol f]~3fl . If the 
evaluation to this accuracy is accomplished, one can judge whether the overcharging condition (0 is satisfied or not. At the same 
time, we have to examine how the absorption condition (O is modified once the back reaction effects are taken into account. 

To achieve this goal, we need to calculate both the particle motion with the self-force and the self-energy of the system to 
0(e 2 ). However, these quantities have never been computed to such a high order even in the case of Schwarzschild background. 
Here, we propose to bypass these difficult tasks by focusing on the marginal orbits that pass through an unstable equilibrium 
state. The basic strategy that we use here is to relate the unstable equilibrium state to an exact solution known as the double 
Reissner-Nordstrom static solution fl4ll (See also Ref. lfl5ll ) 2 . Since the latter is an exact solution, it includes all the backreaction 
effects. Therefore, we can read the total energy of the system in the unstable equilibrium state from this solution. 



2 Within the linear perturbation, the global perturbative solution describing the equilibrium configuration of the RN black hole and a charged particle was 
constructed in Ref. jUjl. 
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Then, we compute the energy emitted to infinity through the infall of the particle from the equilibrium state to the black hole 
horizon to the accuracy of 0(e 2 ), which can be readily achieved using the standard black hole perturbation method. To compute 
the energy flux emitted to infinity, back reaction effects on the trajectory of a particle can be neglected since they give only 
higher order corrections to the energy flux. As a result, we obtain a sufficiently accurate estimate for the total energy of the final 
state for the marginally plunging orbits. In the two successive sections, we prove that the overcharging is impossible for any of 
such marginal orbits. 

To extend our argument to more general cases, we assume that the other plunge orbits result in final states that have higher 
energies than the final state in the case of the marginal trajectory, whose energy is deduced under the conditions that the mass 
and charge of the initial black hole and the charge of the plunging particle are fixed. Once we accept this rather reasonable 
assumption, the avoidance of overcharging in the marginal cases is extended to general plunge orbits. Although, strictly speaking, 
we cannot verify this last statement, we can at least claim that there is no evidence that supports the possibility of violating the 
weak cosmic censorship conjecture in the present context. 



III. TOTAL ENERGY OF THE SYSTEM WITH A CHARGED BLACK HOLE AND A CHARGED PARTICLE IN 

EQUILIBRIUM 

As a first step to examine the overcharging condition (O for the marginal orbits, we here extract the total energy of the system 
composed of a charged black hole and a charged particle in the equilibrium state, E eq , using the double Reissner-Nordstrom 
(DRN) solution of the Einstein-Maxwell equations lfl4l [l5tl . On the other hand, the total charge is simply given by the sum of 
the charges assingned to respective objects as Qtotai = Q + Q since the electric charge is not carried by a field. In this section, 
we will show that the energy at the equilibrium state is always larger than the total charge, E cq — (q + Q) = 0(e 2 ) > 0. 

A. The double Reissner-Nordstrom solution 



We use the form of the DRN solution presented in Ref. 1U4H in the same notation, except for the signature of the metric, 
(— , +, +, +). The DRN solution is characterized by five parameters: masses of the RN sources mi, m 2 , their charges e l5 e 2 , 
and the separation between them I. It should be emphasized that the mass parameters mi and m 2 do not refer to the rest masses 
of the RN sources ifTill . As we explain later, there is one relation among these five parameters, and hence four of them are 
independent. With these parameters, the metric and the vector potential are written in the cylindrical Weyl coordinates as 

o 2 

ds 2 = -H(p, z)dt 2 + F{p, z)(dp 2 + dz 2 ) + —f — -d(j) 2 , (10) 

H(p, z) 

At = -$(p, z), A P = A Z = A^ = Q. (11) 

To express the functions H, F and $ in a simple form, we introduce two sets of bipolar coordinates (n, 0\) and (r 2 , 9 2 ) defined 
by 



p = y/(n - mi) 2 - a'{ sin6>i , \p = \J (r 2 - m 2 ) 2 - a\ sin6> 2 , 
z = z\ + (n — mi) cos#i , I z = z 2 + (r 2 - m 2 ) cos6> 2 , 



(12) 



where z\ and z 2 are set to satisfy I = z 2 — z\ > 0, and af and a 2 are defined by 

o\ = mf — e\ + 2ei7 , er 2 = m 2 — e 2 — 2e 2 7 . (13) 
Using these coordinates, the functions H, F and $ are written as 

H(p, z) = [(n - mi) 2 - a 2 + 7 2 sin 2 2 ] [(r 2 - m 2 f - a 2 + 7 2 sin 2 1 ] D -2 , (14) 

F(p, z) = [(n - mi) 2 - a 2 cos 2 6 X ] _1 [(r 2 - m 2 ) 2 - a 2 cos 2 9 2 ] D 2 , (15) 

<&(p, z) = [(ei - 7)(r 2 - m 2 ) + (e 2 + j)(n - mi) + j(miCos9i + m 2 cos6» 2 )] D^ 1 , (16) 



where 



and the parameter 7 is defined by 



33 = rir 2 - (ei - 7 - 7 cos6' 2 )(e 2 +7-7 cos^i) , (17) 
7 = (m 2 ei - mie 2 )(£ + mi + m 2 y x . (18) 



5 



The parameter a\ is chosen to be negative while erf positive. The object labeled with 1 is a naked singularity that has a critical 
spheroid at {n = mi, < p < \ai\, z = zi}, while the object labeled with 2 is a black hole surrounded by an event horizon 
at {p = 0, Z2 — o"2 < z < Z2 + 02}. The position of the naked singularity is in the region r\ < mi, which is not covered by the 
original coordinates (p, z). To guarantee the existence of an equilibrium state without any conical singularities between the two 
objects, a balance condition, 

mim 2 = (ei - 7)(e 2 + 7) , (19) 

must be satisfied. This gives one relation among five parameters (mi, 777,2, ei, &2, £), which can be satisfied only when o\a\ < 0. 
Namely, the DRN solution should consist either of a black hole and a naked singularity, or of two extremal black holes. The 
formar case requires that o\ and <j\ should have opposite signs, as we have chosen above. The latter is a special case of the 
Majumdar-Papapetrou solution fl7il . 

The total mass and charge of the DRN solution can be read from the asymptotic forms of the metric functions ( fT~4T > and < fT~5b . 
and the electric potential dT6l >. At the space-like infinity, defined by the limit p — > 00 or z — > 00, the two bipolar coordinates 
(?"i,0i) and (7-2, #2) coincide with the spherical coordinates (r,6) defined by p = rs'mO and z = rcos9. Thus, the metric 
functions and the electric potential are expanded as: 

H = 1 _mi + m 2+o n 



„ 1 nil + ?7i2 „ / 1 
F = 1 — — - + 1 



2 



r \r* 



<£ = ^-^ + 0[-), (20) 



in this region. For asymptotically fiat spacetime, the coefficients of the 0(1 /r) terms of the metric function and the electric 
potential correspond to the total energy and charge of the system, respectively. Thus we find that the total energy is mi + m.2 
and the total charge is ei + e 2 . 



B. Mass dominance of the DRN solution 



In this subsection, we study the DRN solution consisting of a black hole with charge 62 and mass m 2 satisfying a\ > 0, and 
a naked singularity with charge e\ and mass mi satisfying a\ < 0. We assume 

?Ti2 > mi > , e 2 > a > , (21) 

which means that the mass and the charge of the black hole are larger than those of the naked singularity. Under these conditions, 
we show that the total mass is larger than the total charge, i.e. mi + m 2 > ei + ei. 
To this end, we re-express the balance equation ( TT9b as 

7(7) := 7 2 + (e 2 - ei)7 + mim 2 - eie 2 = . (22) 

In addition, from a\ < 0, u\ > 0, and using Eq. ( fT3b . we obtain the conditions for the parameter 7 as 

7 < 1 , 7 < ~\ , (23) 

2ei 2e 2 

and, from the positivity of the separation i > and using Eq. ( fTSl l, we find 

m 2 ei - mie 2 



mi + m 2 



< 7 < 0, for TO2ei — mie2 < 0, 



(24) 



?n 2 ei - mie 2 

U < 7 < , tor 77i9ei — mie9 > L). 

777i + m 2 



In this context, the balance condition d22l is interpreted as the condition that the quadratic equation of 7 has a solution in the 
range specified by the conditions, (l23T l and (f24l . 

Now we examine if a solution exists for three cases: (i) 7 < and mim 2 > eie2, (ii) 7 < and m\m2 < e\e%, and (iii) 
7 > and m\m2 < e\e 2 . For 7 > 0, the case with m\mi > e\e2 is immediately excluded by the condition ( f22l . 
Case (i): In this case, the conditions ( f23l ) and ( f24b imply 

ml - el el - ml m x e 2 - m 2 ei 

- 7 > — 7, > -7 > — 7, » -7 < ; ■ (25) 

zei 2e2 mi + m 2 
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The first two inequalities in d25l l lead to 



7 > 



1 



mi + 77l 2 



m 2 (m\ — e\) mi(e 2 — m 2 ) 



2ei 



2e 2 



> 



mie 2 — m 2 ei 

2(m 1 + m 2 ) 
mie 2 — m. 2 ei 

mi + m 2 



eie 2 



where we used mim 2 > eie 2 in the last line. However, the inequality 
Thus, case (i) is excluded. 

Case (ii): In this case, the balance condition d22l i implies 



(26) 



is in contradiction with the last condition in ( |251 ). 



7 2 + (e 2 - ei)7 = eie 2 - m 1 m 2 > 0, 

which leads to the inequality 

7 + (e2 - ei) < . 

Combining this inequality with Eqs. ( [24-b . we find 

mie 2 - m 2 ei 



mi + m 2 



> — 7 > e 2 — ei 



and hence we have 



m-ye-i > m 2 e 2 



(27) 



(28) 



(29) 



(30) 



However, the last inequality miei > m 2 e2 is in contradiction with the assumption m 2 > mi > and e 2 > e\ > 0. Thus, case 
(ii) is also excluded. 

Case (iii): Recalling that 1(0) = m\m 2 — eie2 < 0, the necessary condition that 1(7) = has a positive root satisfying the 
condition (l24l i is 



/ m 2 ei - mie 2 
V mi + m 2 



mim 2 



(mi + m 2 y 



[(mi +m 2 ) 2 - (ei + e 2 ) 2 ] > 



This inequality leads to 



mi + m 2 > ei + e 2 



(31) 



(32) 



Since both case (i) and case (ii) are not allowed, the only possibility is case (iii), and we also find that 7 cannot be negative 
simultaneously. In case (iii) the inequality d32t was proven. Hence, to conclude, under the assumptions m 2 > mi > and 
e 2 > ei > 0, we have shown that 



E cq > q + Q, 



(33) 



which means that the total energy never goes below the extremal bound in the equilibrium configuration. 



C. Remarks on the use of the DRN solution 



In the preceding subsection, we used the DRN solution as such that describes an equilibrium state of a charged particle in 
a black hole spacetime. However, strictly speaking, we also have to prove that the charged singularity in the DRN solution 
is a valid approximation for a charged particle. An arbitrary solution that allows a singularity may in general possess higher 
multipole moments other than the monopole. In that case the charged particle will have an extra self-energy due to higher 
multipole distortion. Thus, one may suspect that there might be another equilibrium configuration that has lower total energy 
with the same total charge. Hence, it would be necessary to explicitly demonstrate that the charged particle in the DRN solution 
does not possess higher order multipole moments that are large enough to contribute to the total energy of 0(e 2 ). A detailed 
discussion is given in Appendix lAl and the outline is as follows. We first show that the behavior near the singularity in the 
DRN solution can be described by a perturbed RN configuration. Perturbative expansion does not break down even near the 
singularity in the DRN solution. From this expression we can read the maximum amplitude of the multipole moments that the 
particle possesses. The largest component is the dipole and its magnitude is 0(e 3 ), and we find that the contribution of such a 
dipole moment to the total energy is at most 0(e 3 ). 
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IV. ENERGY FLUX RADIATED FROM A CHARGED PARTICLE 

In the preceding section, with the help of the DRN solution, we showed that the total energy of the equilibrium state of a 
charged particle in a RN black hole cannot be less or equal to its total charge. This suggests that the system in equilibrium 
does not exceed the extremal bound, but it is not sufficient yet to exclude the possibility that a RN black hole is overcharged 
by absorbing a charged particle. The effects of the electromagnetic and gravitational radiation can reduce the total energy as 
the particle falls from the equilibrium position to the event horizon of the black hole. Then the total energy in the final state 
will be given by £" cq — E x with E x being the energy radiated away to infinity. In this section, we estimate using linear 
perturbation theory. We examine if the signature of E cq — q — Q — E^ can be negative, to see whether the overcharging in the 
final configuration is possible or not. 

A. Correspondence between two different pictures 

Before evaluating the radiation energy emitted to infinity, we would like to briefly mention the correspondence between the 
DRN solution and the equilibrium configuration of a test particle in a RN black hole. Here we refer to the latter as the black 
hole perturbation (BHP) picture. Our current interest is in the case in which the charge of the test particle is 0(e), while the 
difference between the total energy and the total charge is 0(e 2 ): 

e x :=ae, e 2 := 1 , m 1 :=0{e), m 2 := 1 + zu(l + 0(e)) , (mi + m 2 ) - (e x + e 2 ) = 0(e 2 ) . (34) 

As a starting point we do not exclude the possibility that w is as large as 0(e). There might be more subtle cases if we consider 
the possibility that the difference (mi + m 2 ) — (e\ + e 2 ) is 0(e 3 ) or higher, but such cases are beyond the scope of the present 
paper. 

Then, as shown in Appendix|B] the parameters used to describe the DRN solution are expanded in terms of e as 

ei = ae , e 2 = 1 , 7 = .ge 2 , £y = ge 3 , (35) 

where a and g are expected to be 0(1), while g can be either 0(1) or even smaller. If g is of higher order in e, I should be large 
so as to maintain £j = 0(e 3 ). The meaning of a is obviously identical to that in the BHP picture. Substituting Eqs. (l35T l into 
Eqs. (fT3T l, ([T8l and ( fT9] l, we obtain approximations for the mass parameters (mi, m 2 ) and (01, er 2 ) as 

2 $ 3 1 i / i 9 \ 2 2(2 ~\ 4 2 S 2 /o/r\ 

mi = ae-ge - -e , m 2 = 1 + I g + — I e , a x = (g - ag)e , <r 2 = -e , (36) 

where the terms of higher order in e are neglected. The relations ( |36T l enable us to evaluate the difference between the total 
energy and the total charge as 

7m+ m 2 -(ei+e 2 ) = |-e 2 + 0(e 3 ), (37) 
la 

which is positive as long as 7 and I are positive. This is, of course, consistent with the more general argument given in the 
preceding section. 

The relation between the parameters g and g in the DRN solution and the parameters in the BHP picture is examined in 
Appendix ICl and is given by 

g = 4a + 0(e), g = 2\/a 2 - c 2 + 0(e) = 26 + 0(e). (38) 

(Here b is simply determined by the equilibrium condition (ICU .) From these expressions we can confirm that there are counter- 
parts of g and g for any values of a and c. Using the above relation, we finally obtain 

7771+ ?77 2 -(ei+e 2 ) = 2e 2 + 0(e 3 ). (39) 

In the analysis neglecting all back reaction effects, as in Ref. 0], our interest was restricted to the parameter region with b > 1 
that allows overcharging in the final state. Under this restriction, g must be definitely O(l). As a result, I must be O(e). As 
shown in Appendix. ICl the corresponding proper separation is O(l). However, now we know that, once all the back reaction 
effects are taken into account, the difference between the total energy and the total charge at the equilibrium configuration is 
always positive and of 0(e 2 ) in our setup of the problem. Hence, if the effect of energy loss due to radiation is less than 
0(e 2 ), there is no possibility of overcharging. To the contrary, if the energy loss is as large as 0(e 2 ), all parameter region in 
principle may give an example of overcharging. Hence, in the succeeding discussion we will not restrict ourselves to the case 
with b = 0(1). From the above correspondence, we understand that, when b is of higher order, so is g, which corresponds to 
the case that I is greater than O(e). In this case, the difference in the rest mass and the charge of the point particle a — c is also 
higher order and the proper separation between the horizon and the particle is as large as 0(| log e|). 



B. Gauge invariant perturbation of the RN black hole excited by a charged particle 



To estimate the energy emitted to infinity by the electromagnetic and gravitational radiation, we use linear perturbation theory, 
which has been developed for a RN spacetime in Refs. lfl2l [l8ll . Here we adopt the formulation developed by Kodama and 
Ishibashi lfl2ll . in which the equations are reduced to a set of decoupled second-order differential equations. 

Owing to the background spherical symmetry, one can decompose the perturbations into scalar- and vector-types in general. 
It is easy to check that the energy-momentum tensor of a charged particle in a radial orbit is purely of scalar-type, and therefore 
only the scalar-type perturbations are excited. In addition to that, by choosing the coordinates so that the trajectory is along 
the axis we can restrict our analysis to the axisymmetric mode m = 0, where m is the azimuthal index of spherical harmonics 
Yi m (9, <j>). As shown in Ref. IU2I1 . the equations of the scalar-type perturbation are reduced to two decoupled equations for two 
gauge invariant variables, $±. They are given in terms of Fourier-harmonics expansion as 



dr* 2 



+ (uj 2 -V ± (r))<Z> ± (r*) = S ± (r*,u), 



(40) 



where the tortoise coordinate r* is defined by dr^/dr = l//(r), when f(r) is as defined in Eq. Q, and us is the frequency. 
For simplicity, we abbreviate the labels [, m, u to be attached to the master variables, effective potentials and source terms. The 
effective potentials, V± (r), are defined by 



V±(r) :-- 



f(r)U±(r) 
64r 2 ff 2 » ' 



with auxiliary functions 

H+ := 1 — 3Sx , H- := to + 3px , S 



2m \M ) ' 9 



1 



p := 1 + mS 



(41) 



(42) 



where x := 2M/r and z := Q 2 /r 2 . Here we introduced m := [(I + 1) — 2 where [(I + 1) is the eigenvalue of the Laplace 
operator on S 2 . (It should be stressed that m is different from the azimuthal index m). The functions U± (r) in Eq. fiTt are given 
by 



U+{r) := -25925V 4 + 576(3to(5 + A)S 2 x 3 - I926x{3x + m) + 64(m + 2) , 
U_(r) := -25925p 3 x 4: - 576(3to<5 - l)p 2 x 3 + I92mpx(3x + m) + 64m 2 (to + 2) 



(43) 



The tensor composed of the scalar-type harmonics vanishes for [ = 1, or equivalently for m = 0. This mode is called an 
exceptional mode in Ref. 03], which requires special treatment. In this case, although the master equation d40b is still valid, 
the dynamical variable is <£> + only. For the exceptional mode, 5 takes the value Q 2 /9AI 2 , which is obtained by taking the limit 

TO — > 0. 

The source terms S± (r* ,u>) are constructed from the energy-momentum tensor and the current induced by a charged particle. 
They are given by 



S±(r*,w) := a±(r)5*(r,w) + b±S A (r,w), 



(44) 



with the coefficients defined by 



K(r), 6+} 
{a_(r), &_} 



mQ 3(M + u)Q 3(M + u)k 
~2~ + 



2r 



V2 



3{M + v) 



-4\/2Qk| , 



(45) 



and 



S*(r, lo) 



V2fQn 



Psi 
H 



rH \ H 



t_ 

iiO 



2(to + 2) 

i dsr 



iui 



2r 



Or 



Sa(t, us) := - 



r 2 H 



^--fi-UJr)- 

ius or lustiK r 



2r 2 S r r 
2y/2fQ S 7 t 



(46) 
(47) 
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Here functions H, P$i and Pg 3 are defined by 

H :=m + 3x-iz, P S x := 4x[2z - 3x + 6 + m(m + 4)] , P S3 := 4{3x - 8z) , (48) 

Sa(t) and J a (i) are respectively the Fourier components of the energy momentum tensor and the charge current of a radially 
falling charged particle that starts with a stationary point r = ro in the infinite past t — > — oo. Their explicit forms are expressed 

as 



ST = -«W I —Yio 0(r - r), S r r = ^/x/" 1 ^ ) 



Jt = , , e lul F ra 0(r o - r). J r = — — ^— — e wl Y [0 9(r - r), (49) 
[(1+1) v ; ' [([+1)V^/ 

where Yio = ^/ (21 + l)/47r. In these expressions dT/ds, dR/ds, dT/dR :— (dT / ds) / (dR/ ds) and T are to be understood as 
functions of R, and R is replaced with r. The function T is obtained by integrating Eqs. © and (O and eliminating s. 

The formal solutions <1>± of Eq. d40b applicable to both the generic and exceptional modes are obtained by the usual Green 
function method. The asymptotic form of the solution at r — > +oo is given by 

*^">=TOtr (50) 

with 

X±(w) := f ° $5(r*,a;)5 ± (r*,u;)dr*, (51) 



where the functions $£ P ^ ln (r*,u;) are the homogeneous solutions of Eq. (|40b that satisfy the boundary conditions $± p (r*, uj) — > 
e iur at r * oo an( j $ 1 _g(7-* ; w ) — >• e~ lur at r* — > — oo, respectively. We also introduced the Wronskian defined by 
W[$\$ 2 ] := (<9 r .$>*))$ 2 (r*)-$V*)(dr*$ 2 (r*)). 

The energy flux carried by electromagnetic and gravitational waves to infinity is described in terms of the variables $± 
as given in Eq. ( ID5I ) in Appendix iDl In this expression $±(r, t) are the complex conjugations of $±(r, t). Substituting the 
solution ( T5Tb into the time integral of Eq. dP5t , with the aid of Parseval's theorem, we obtain the total energy radiated to infinity 



87rl([+l) w 2 f\X + \ 2 , (I-l)([+2) |X_| 2 , 



v{M + v) \\W + \ 2 16 \W-\ 



where W± := VK[$^ P , The explicit form of X± is given in Eq. dP9b . We point out that the integrals in the expressions for 
X± given in Eq. dP9t , with the last two terms replaced with Eq. (1D1 U all take the form 

I(r)&£(r)e iuT( - r Ur, (53) 

where 7(r) is a certain regular function of r, whose typical scale of variation Ar is 0(M~ 1 ). For large tu, the integrals d53l take 
the structure of rapidly oscillating function ^ 1 ^(r)e lwT ^ r ' multiplied by the slowly varying function I(r). In general, such an 
integral is known to be exponentially suppressed for large Mu>. Thus, the w-integral in Eq. ( |52l has an effective high frequency 
cutoff at 0{M~ X ). 

Now we show that this radiated energy is 0(e 4 ) or higher. It is obvious that |X±(w)| is suppressed by a factor of O(e) since 
each term manifestly contains \i or q. Therefore what we have to show is that there is an additional suppression of 0(e) in 

|x±H|. 

As a preparation, we show that the inverse of the Wronskian is suppressed in the limit oj — >• like W± 1 oc (Mw)'. The 
presence of this suppression is understood as follows. Since the Wronskian is constant independently of r, we evaluate it in the 
limit r — ► oo as 

iuj&?{r) e 2ur . (54) 

or J 

Note that in the limit r — > oo, the difference between r* and r is suppressed by 1/r, so that we can use r instead of r* in Eq. d54b . 
When r is large enough compared with M, the homogeneous equation corresponding to Eq. d4Qb can be approximated as 
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The general solution of this equation, which is written in terms of the Bessel function of the first kind J, and that of the second 
kind Yu should also describe $+. Hence, we have 



$± * CjJ^fj l+1/2 (nj) + C Y J^Y l+1/2 (ru>) , (56) 

where Cj and Cy are coefficients to be determined by the condition imposed near the horizon. Then, for a small r of 0(M) 
on the verge of the validity of this approximate solution, the two terms in Eq. d56*b should be equally important in general. 
Furthermore, as there is no significant feature in the potential, the amplitude of these terms must be O(l). This determines the 
order of magnitude of the coefficients as Cj = O^AIuj)" 1 ^ 1 ) and Cy = 0((Mlu) 1 ), since for a small roj the above expression 
asymptotically behaves as 

7r / cor \ 1+1 1 1 fujr\- 1 . 



^-^/aU-J WTW)- Cy V^^ ] r( - [ - 1/2) - (57) 



On the other hand, for a large r, the first term dominates to give 



7T 



&l -> Cj cos [ur- -(l+l)) . (58) 

Substituting this expression into Eq. (1541) . the order of magnitude of the Wronskian is estimated as W± = O (j(M(j)"'). As 
anticipated, we find that the inverse of the Wronskian scales like oc (Mw) 1 . 

Now, in order to prove the presence of an additional suppression factor in |X±(w)|, we focus on the fact that (dR/ds) 2 = V^-°' 
is always suppressed in the present setup. It will be obvious that V^°'(r) is a quadratic function of 1/r bounded from below. 
Furthermore, we know that both V^°' and dV^ jdr vanishes at r = Tq. Therefore takes its maximum value at r = r + , in 
the interval of our interest between r + and r$. Hence, we have 

V^\r)<V^{r + ) = \(E qQ ) = l^_f}! e 2 + (e 3 ). (59) 

A' 2 V M + VM 2 -Q 2 J c 2 

This implies that the velocity of a particle dR/ds given in Eq. (0]i is always at most 0(e). When the particle moves very slowly, 
the amount of emitted radiation is also expected to be small. This intuition can be made explicit in the expression for X± as 
follows. One can see that each expression in Eq. d53l ) is regular on the boundaries of the integral. Here we replace e lwT ^ with 
an equivalent expression 

and perform integration by parts. The regularity on the boundaries of the integral is not ruined thanks to the presence of a factor 
/ and dR/ds in the above expression d60l ). (Notice that the combination f(dT/ds) is regular on the horizon.) This manipulation 
adds at least one dR/ds factor at the expense of decreasing the power of u> by one. Repeated application of this integration 
by parts is restricted by the requirement for the convergence of the w-integral in Eq.d52ll. However, this is not such a severe 
constraint, because the inverse Wronskian squared IW^I -2 gives suppression for a small u> proportional to uj 21 as mentioned 
above. Therefore this additional suppression owing to the dR/ds factor guarantees that the total energy emitted to infinity is at 
mostO(e 4 ). 

To summarize, the total energy emitted to infinity by a particle that falls from the unstable stationary point to the horizon, Eqo, 
is always suppressed, and it cannot be as large as 0(e 2 ) for any parameter choice. Hence, the effect of the energy loss due to 
radiation does not affect the inequality (|37l l, which is a relation at the level of 0(e 2 ). Thus, we conclude that 

E cq - Eoo > q + Q, (61) 

is always satisfied. Namely, the total mass of the final state can never be reduced below the critical value that corresponds to the 
extremal bound. 



V. CONCLUSION 



In this work, we have examined the back reaction effects of 0(e 2 ) when a charged particle whose mass and charge are of 
0(e) is absorbed by a nearly extremal RN black hole. To avoid the technical difficulties related to the electromagnetic and 
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gravitational self-force, we concentrated on the case of the marginal orbit, which is the separatrix dividing the plunge and recoil 
orbits. We first showed that, with the aid of an exact solution, the total energy of the system is always greater than the total charge 
for the equilibrium configuration that the marginal orbit passes through. Then, we demonstrated that the radiative energy loss as 
the particle is falling into the black hole from the equilibrium position is 0(e 4 ) or higher. Combining these results, we succeeded 
in proving that the total energy of the system composed of a charged particle and a black hole is always greater than their total 
charge for the marginal orbit, once back reaction effects are properly taken into account. In short, the back reaction effects 
prevent a nearly extremal RN black hole from being overcharged, and hence we conclude the cosmic censorship conjecture is 
not violated. 

As was mentioned in Sec. MB I our discussion relies on the assumption that the marginal orbit passes through the unstable 
stationary configuration. Although it is difficult to imagine that this is not the case, our modest conclusion at the moment will be 
that there is no evidence that supports the possibility of overcharging a nearly extremal RN black hole by absorbing a charged 
particle. If we further accept a rather natural assumption that the final state of all the plunging orbits has a higher energy than 
the case of the marginal orbit for given any charges of particle and black hole, we can exclude the possibility of overcharging 
without the restriction to the marginal orbit. 

To obtain a definitive answer to the question whether these assumptions are correct or not, it is necessary to solve the orbital 
evolution directly including the self-force effects. The formal expressions for the gravitational self-force in vacuum spacetime 
was already derived by Mino, Sasaki and Tanaka, and independently by Quinn and Wald [ 1911 . and is known as the MiSaTaQuWa 
force. In the case of a charged particle in the non-vanishing electro-magnetic background, even the fundamental formulation 
corresponding to the MiSaTaQuWa force is lacking. Although there is a recent work in this direction I2OII . it has not yet been 
developed to the level applicable to the overcharging problem discussed in this paper. 

The main lesson of this paper is that it is important to include the non-dissipative part of the self-force when we examine 
the possibility of overcharging. In the context of the spinning up of a Kerr black hole proposed in Ref. @J, Barausse et al. 
have shown that, to avoid overspinning, it is not sufficient to take into account only the dissipative part of the self-force, i.e. the 
energy loss and the angular momentum loss due to gravitational radiation ifTTIl . As was discussed in Ref. Hill , also in this case 
it is essential to take into account the non-dissipative part of the self-force. However, in contrast to the overcharging process 
discussed in this paper, we cannot expect the existence of the stationary intermediate configuration characterizing the marginal 
orbit. Therefore, directly analysis of the self-force will be unavoidable. 
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Appendix A: Geometry near the singularity in the DRN solution 

In Sec. lHlCl we implicitly assumed that the singularity in the DRN solution is the counterpart of a charged particle. The aim 
of this appendix is to confirm this correspondence. The main focus is on the asymptotic behavior near the singularity, so as to 
show that the deformation from the spherically symmetric singularity is sufficiently small. 

To analyze the geometry near the singularity in the DRN solution, the original (p, z) coordinates are inappropriate since they 
do not cover the region where the singularity resides. In order to analytically continue DRN metric ( [10) . we can use the bipolar 
coordinates (ri,#i) defined by Eqs. (112) . However, to see the behavior near the singularity, it is more convenient to change 
the radial coordinate from ri to D. Furthermore, as we are focusing on the behavior near the singularity, we use the rescaled 
variable A defined by D — Ae 2 , instead of V. For simplicity, we concentrate on the case that the stationary point is rather close 
to the black hole horizon, in which g is 0(1). In the case with a smaller g, which corresponds to a larger I, the deformation is 
expected to be even weaker. 

When g is 0(1), it is more convenient to use the parametrization A = g/g instead of g introduced in Eqs. (135) . Hence, our 
parametrization adopted here is 

ei = ae , e2 = 1 , 7 = ge 2 , t = Ae. (Al) 

The other parameters mi,m2, a\ and a\ are to be determined from Eqs. (113) . dT~8b and < fT9b to the sufficiently high order in e. 
For the present purpose, the approximate expressions ( f36t are not sufficient. 
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To write down the metric and the electric potential in the (A, X ) coordinates, we need to describe r±, r 2 and 2 in terms of 
(A, 9\). From the definitions (ITZb and ( fTTI i, setting z\ = and z 2 = I, we obtain 

[r\ — m 1 )cos0 1 — I 
r 2 = m 2 H . 



sin 2 2 



ri 



cos 02 

[(ri - mi) 2 - cr 2 l sin 2 0i 



(r 2 - m 2 ) 2 - cr| 
1 [( e i — 7 ~ 7cos02)(e2 + 7 — 7cos0i) 



(A2) 



which can be solved by iteration in the presented order, once we have appropriate initial values for r\ and 02. The location of the 
singularity is specified by V — 0, where the lapse function g tt = —H and the electric potential $ diverge. This point is expected 
to be close to the critical spheroid, and hence we have n ~ mi. Then, we also have 02 ~ — ir near the singularity. Using these 
crude estimates as the initial condition for the iteration, we can solve these equations to a sufficiently high order in e. 

Since the coordinate A is not best suited for seeing how the DRN solution deviates from the single RN singularity, it is better 
to further replace the radial coordinate to a physically well-motivated one. One simple possibility is to adopt the inverse of the 
electric potential as the radial coordinate, i.e. 



R := 



act 



where at is the charge of the particle and we will find that the factor ct takes care of the change of the time coordinate, 

dT := ctdt, 



(A3) 



(A4) 



with c being a constant of 0(1) to be determined later, once we notice that the combination A t dt(= —<&dt) is gauge invariant. 
Here, the new time coordinate T appropriate to describe the charged particle as a perturbed RN geometry elapses slower than the 
time of an asymptotic observer t by the red shift factor ct since the particle is staying deep inside the gravitational potential of 
the black hole. Also, a constant $ was introduced to adjust the zero point of the electric potential. $ w iH be also determined 
later. To focus on the range where c\jR = at/ R = 0(1), we mainly use the rescaled radial coordinate g = R/t, assuming 
g = 0(1). (Notice that the coordinate g is different from g in the cylindrical Weyl coordinates, used to describe the DRN 
solution in Eqs. (TTOb and (fTTT i.) 

Then, it is straightforward to expand the metric in these coordinates up to O(t). By choosing 



y/a\(a\ - g) X(g 2 - 2ag\ + 2a 2 A 2 ) + 2 



$ = 1 - Xt 



the resulting metric is expressed in the following form 



a 2a y / aX(aX - g) 

2 X(2aX-g) t ^ 3 



2a 



+ 0(6 J ), 



(A5) 



9TT = -f(Q) 

1 

9RR - JTT 
9ij = r2 Hj 



1 — t COS 01 



1 + t COS ( 



- +0(t 2 ) 
\^aX(aX-g)J 

6(a 2 + g 2 ) 2g(6aX-5g) 

a y aX(aX — g) 



+ 0(e 2 ) 



i + 2eC os0i { A(a2 ± g2) - W^A^) + (t 2 ) 



9R8! = eR 2 sin 0i 



aX 



a A 

2 > 



a V {aX - g) 



+ 0(£ 2 ), 



(A6) 



where we have defined 



2 v / Q A(aA-, 9 ) ^ g{2aX-g) \ 1 1 = f_ _ 2ji 1 
A 6 v /aA(aA - g) J Q R 2 R 



(A7) 



The second equality gives the expressions of q and fj, in terms of a, g, X and e, which are consistent with q = at, \i = ct in 
Eqs. (O supplemented with the relations Xg = g = 4a + 0(t) and g = 2\J a 2 — c 2 + 0(t) in Eqs. d38l l. In the next leading 
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order in e the metric given in Eqs. JA61 > has only dipole type perturbation at around R = 0(e). Perturbations belonging to higher 
multipoles are of 0(e 2 ) or higher. 

The asymptotic behavior of perturbations for small R is to be determined in line with the nature of a charged particle. For 
example, when an external dipole field is imposed, a charged particle will gain some induced dipole moment. However, the 
magnitude of the induced dipole moment will depend on the property of the particle to some extent. In this sense there remains 
arbitrariness in the most relevant choice of the boundary conditions for our charged particle. Instead of discussing this, we show 
that the effect of the induced dipole of 0(e) in the expressions ( IA6b does not contribute to the total energy of the system at the 
level of our interest, 0(e 2 ). 

The key observation is that it is at R = 0(e) where the dipole perturbation is O(e). This dipole field is a superposition 
of the external field caused by the presence of a charged black hole and the induced dipole of the point particle. Hence, the 
latter amplitude will be also at most 0(e). The dimensionless perturbation caused by an object having a dipole moment D is 
0(D I R 2 ). Therefore the magnitude of the dipole moment is evaluated as D = 0(e 3 ). In the region R ^> 0(e), the configuration 
can be understood as a RN black hole solution perturbed by a charged particle and an electric dipole placed at the same place. 
In such a regime the effect due to the dipole is at most 0(e 3 ). On the other hand, in the vicinity of R = we can evaluate 
the energy due to the dipole from the expanded metric presented in Eqs. ( IA6b . The explicit expanded metric shows that the 
monopole part of the perturbation starts with 0(e 2 ). Since the change of the mass 6fj, will appear in the monopole part of the 
dimensionless metric perturbation as 6/J./R, we find that 5(1 is at most 0(e 3 ) in the region R < e. As a result, we conclude that 
the contribution of the dipole to the total energy is 0(e 3 ). The higher multipoles can be discussed in a similar way, and we find 
that they are even higher order in e. 



Appendix B: e-expansion of the parameters in the DRN solution 

This appendix is devoted to deriving the appropriate order of 7 and £ in terms of e-parametrization. Our starting point is the 
assignment of the ordering that defines the setup of our problem given in Eqs. ( f34b . The parameters are also constrained by the 
inequalities a\ < and a\ > 0, which are required because a 2 is positive for a BH and negative for a singularity. 

From the ordering given in Eqs. d34b , the leading term of 7 in Eq. (TT~8T > can be read as 

7=^^(1 + 0(6)). (Bl) 

Substituting the above approximate 7 into the balance condition ([T9T i, we find 

£(e 1 -m 1 ) = 0(e 2 ). (B2) 

Multiplying £ on both sides of Eq. ( IBU . with the aid of (IB2t we have £7 = 0(e 2 ). Next, using the last equality in Eqs. J34t . we 
find 

ex - mi = w + 0(e 2 ). (B3) 

Plugging Eq. (|B3t into Eq. (B2i . we also have £w = 0(e 2 ). 

We now show that Eq. (IB2b is in fact further suppressed. Substituting Eq. (fT8T l into Eq. (TT~9b and expanding it up to 0(e 2 ) 
after multiplication of the factor (£ + mi + 7712) on both sides, and using Eq. (IB2b and the relation £w = 0(e 2 ), we obtain 
2mi(ei — nil — to) + £(&i — mi) = 0(e 3 ). With Eq.(|B3l, this reduces to 

£(ei - mi) = 0(e 3 ), (B4) 

which also implies £7 = 0(e 3 ) in the same manner as above and £w = 0(max(e 3 , £e 2 )) with the aid of Eq. ( |B3t . From Eq. ( fT8] l, 

(l + ro)(ei-7)-l7 

mi = — , (B5) 

1 + 7 

follows. Substituting this expression for mi into erf given in Eqs. dT3l ). with the aid of the relations (IB lb and dB3b . we obtain 

a\ = 7 2 + 0(e 4 ). (B6) 

Thus, the requirement that a\ < can be consistent only when 7 = 0(e 2 ). Then, immediately both ei — mi and zu turn out to 

beO(e 2 ) 

To summarize, if £ is 0(e), we find that both ei — mi = 7 + 0(e 3 ) and w are 0(e 2 ). When £ is not as small as O(e), the 
order of ei — mi and 7 depends on the magnitude of £ in such a way that £(ei — mi) = £(£ + 1)7 + 0(e 4 ) becomes 0(e 3 ). By 
contrast, w stays 0(e 2 ) irrespective of the magnitude of £. 
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Appendix C: Mapping an exact static solution to an equilibrium configuration of a charged particle in the RN spacetime 



In this appendix we consider the mapping of the DRN static solution to an equilibrium configuration of a charged particle 
in a RN black hole background: the BHP picture. It is not so trivial to find the relations between the different descriptions, 
especially the relations between g and g in Eqs. fl38l l, and the parameters a, b and c given in Eqs. (0. To obtain these relations, 
we focus on two geometrical quantities: one is the proper distance L along the symmetric axis between the event horizon of 
the black hole and the equilibrium position of the particle, and the other is the proper distance D from the event horizon to the 
asymptotic region, measured along the symmetric axis in the direction opposite to the particle. For our current purpose, in linear 
perturbation theory analysis, it is enough to estimate the emitted energy at the leading order in e since the energy flux d52l is 
always 0(e 2 ). The back reaction effects on the orbit contribute only to the higher order corrections in the energy flux. Therefore 
we consider the mapping at the level of linear perturbation of a RN black hole, restricted to the order in e expansion necessary 
to specify the setup in the BHP picture. 

We first compute the above geometrical quantities in the BHP picture. Since we consider the equilibrium state, however, the 
parameters a, b and c in Eqs. (0 are not independent but satisfy the relation 

b= Va 2 -c 2 + 0(e), (CI) 

which is derived from the equilibrium conditions dR/ds = dR 2 /ds 2 = 0. The value of the radial coordinate at the equilibrium 
position is also derived from these equilibrium conditions as 

r = 1 + -j^= e + 0(e 2 ). (C2) 
V a c z 

Then, the quantities mentioned above are calculated in the BHP picture as 

r° r a dr' 1 a + c 

L= Vf^dr'= - 7 _ = -log - + 0(e), (C3) 

Jr+ Jr + s/W) 2 a-c 

D(r) = / ^JuTrdr 1 = f -Jl= = - log(e) - 1 - (r + logr) + O (-, e) , (C4) 



vTF) 



where f(r) is the metric function of the RN black hole defined by Eq. (01. 

Next we evaluate L and D in the DRN solutions. In the DRN solution ([T0]>-(fT8]>, we can choose (zi, z 2 ) = (— I, 0) without 
loss of generality. For this choice, the event horizon on the axis is at (p, z) = (0, ±<t 2 )- Then, the proper distance L in the DRN 
solution is formally given by 

P — CT2 

L= y / F(0,z)dz, (C5) 



where F(p, z) is one of the metric functions of the DRN solution defined by Eq. ( fT~5T > and we take the upper end of the integral 
at z = —I, which is sufficiently close to the singularity of the body 1. In the range of — £ < z < —o^ on the symmetric axis, 
from Eqs. (I121 l. we obtain 

\n=z + mi+l, 6»i = 0, 
yr 2 = -z + m 2 , f) 2 = 7r , 

and then the metric function on the axis is given by 

w(n , A _ [{-z + m 2 ){z + m 1 + i)-e 1 e 2 ] 2 
[U,Z> [(z + 1) 2 - a\\[z 2 - a 2 ] ■ K ' 

The direct integration of Eq. ( IC5b gives a complicated combination of incomplete elliptical integrals that is inconvenient for our 
analysis. To avoid the complication, we divide the range of the integral [—£, —a 2 \ into [—£, —£+Ae 3 ^ 2 ] and [— £+Ae 3 ^ 2 , — 02], 
where A is an arbitrary constant of O(l). In the former range, we expand the term 1/ \J z 2 — a 2 in the integrand of Eq. ( IC5b 
with respect to z + £ — (9(e 3//2 )(<C 02), keeping the term 1 / y/(z + £) 2 — u\. On the other hand, in the latter region, we expand 
1 / yj(z + £) 2 — o\ with respect to o\ while the factor 1 / \J z 2 — a 2 is kept unexpanded. Then, one can perform the integrals to 
obtain an approximate estimate of L as 





y/F(0, z)dz = arccosh* / + 0(e) . (C8) 
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The proper distance D from the event horizon of the black hole to the asymptotic region z S> 1 along the symmetric axis is 
also formally given by 

D= lim / y/F(0,z')dz'. (C9) 

*J (72 

Here it should be noticed that, unlike the case of L, we take the lower end of the integral at z = a 2 because the integral is 
performed in the direction opposite to the body 1 (singularity). In the range of the integral on the right hand side of Eq. dC9t , the 
bipolar coordinates are given by 

ir 1 =z + m 1 +£ 1 0i=O, 
\r 2 — z + m 2 , 6 2 = . 

Substituting Eqs. (IClOb into Eq. H5i . the metric function is rewritten as 

F(0 z) = ^ + ^)(^ + m 1+ l)-( ei -2 7 )e 2 ] 2 (z + m 2 f 
( ' ' [{z + if - <jI)[z 2 - a 2 2 ] z*-o\ V + U W>- ^ Ll) 

Then, the integral dC9l ) with Eq. jC lib can be evaluated. Since the coordinate z is related to r on the axis as z = r — 1 + 0(e), 
which can be easily verified by looking at, say, the lapse function, we finally obtain 

D = -]og( ^\~- ) -l + (r + logr)+o(-, eY (C12) 



2 V a I ' \r 

Comparing Eq. (IC4b to Eq. dC12b and Eq. dC3b to Eq. dC8l ). we establish relations between the two sets of parameters presented 
in Eqs. d38l l, and with the aid of the constraint equation (ICU valid up to 0(e) we obtain b. 

Appendix D: Evaluate the energy radiated to infinity from a moving charged particle 

In this appendix we evaluate the energy Eoo emitted to infinity from a charged particle that falls into a RN black hole along a 
radial orbit. 

1. Energy flux formula 

First, we derive the effective energy flux generated from the scalar-type perturbation in the Kodama-Ishibashi (KI) formalism. 
The energy flux carried by gravitational waves and that by electromagnetic waves are decoupled in the asymptotic region r —> 
+00, and they are, respectively, given in terms of the perturbation of the metric h^ iv and that of the electro-magnetic field strength 
/^asiH 

pGw ._ 2 f dQ dh a pdh af - EM 2 f 

h oo —-1 ] J^—Q t <^T> h oo - ~r J dU faff , (Dl) 

in the transverse-traceless gauge defined b y M = 0, V u h^ u = 0. Here, averaging over several wavelengths of radiation 
is assumed. We rewrite these expressions (ID 1 b in terms of the master variables $± in the KI formalism. To do so, as an 
intermediate step we consider $ and A that are related to $-t at the level of coefficients of the spherical harmonics expansion by 

$ ± = a± $ + b±A, (D2) 

with the coefficients defined in Eqs. d45l l. where the indices of spherical harmonics, [ and m, are suppressed for notational 
simplicity. At the leading order in the limit r —> +00, the master variables $ and A are, respectively, related to the perturbations 
of gravitational and electro-magnetic fields as 

hij := [ (l + 1) r^Sij, fai w VW+T)e ab (D b A) , (D3) 

[,m l.m 

where a, 6-indices run over t, r-coordinates, while i, j-indices over angular coordinates. The other components are of higher 
order in 1/r. The totally anti-symmetric symbol e a b is defined so that e tr = 1. The harmonics are defined as 

§ := Y [m , Si := 7 =L=AHm , §ij ■= -T^—D.b.Y^ + ^Yim , (D4) 

y 1(1 + 1) K l + *-) z 
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where Y[ m are the usual spherical harmonics on a unit sphere S 2 and Di represents covariant differentiation with respect to the 
metric of S 2 , %j. 

Then, substituting Eqs. iD2i , (ID3b and dD4b into Eq. ( IDU , the total energy flux escaping to infinity is evaluated as 



E 



l,m v 



E 



4tt[([+1) 



Qk 2 v(M + v) \ dt dr 



8k 2 dt dr 
d<S>+d^>+ ([- l)(t + 2) 
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(D5) 



where quantities with " ~ " represent the complex conjugations, and all terms suppressed by 1 jr are neglected. Note that this 
expression is valid for all modes including the exceptional mode with [ = 1. However, for the exceptional mode the second 
term in the parentheses on the last line vanishes. Hence, we find that $_ is irrelevant for this mode, which is consistent with the 
fact that physical degrees of freedom for the metric perturbation are absent for the exceptional mode. (See Ref. fl2ll for further 
details.) 

We would like to add a short remark about the expression for the energy flux ( ID5b . This is a special case of a more general 
expression for the energy flux defined as an integration over an arbitrary 2-surface of constant r: 



4tt[([+1)/ /d$+d$+ (t- l)(I+2) <9$_ 



9k 2 v(M + v) \ dt dr 
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dt dr 



(D6) 



which is identical to Eq. (lD5b in the limit r -> oo. Here, / is the metric function of the background RN black hole given in 
Eq. ([3]). To verify that the above expression is conserved in the absence of source terms, we focus on the fact that the master 
equation d40b is derived from the variational principle of the action 



(D7) 



where gfy is the two dimensional metric composed of the (t, r) components of the RN metric, vi 2 ^ is the covariant differentia- 
tion with respect to gfy and <j(2) is the determinant of this two dimensional metric. The t and r components are labeled by the 
Latin indices, a, b. Once the action (ID71 > is at hand, its variation with respect to g?K gives a symmetric tensor T a b that satisfies 

the conservation law V? T a 6 = 0. Then, one can construct a conserved current T a b(dt) b associated with the background Killing 
field (dt). Thus, we find the conserved flux per unit coordinate time is proportional to 



T ab n a (d t ) b 



„d$± <9<J> d 



(D8) 



where n a is the unit vector normal to a surface of constant r. Since the expression iD6i is a linear combination of these conserved 
fluxes, it is guaranteed to be conserved. For any wave packet the time integral of this flux is independent of r, and agrees with 
the net energy flux evaluated at infinity. The net energy flux is also a quantity that must be independent of r. Therefore for any 
wave packet the time integral of this flux (ID6t evaluated at any radius gives the net energy flux, and we conclude that Eq. ( ID6b 
is indeed an expression for the conserved energy flux valid for any value of r. 



2. Evaluation of the energy flux from a radially falling particle 

Now we turn to the evaluation of the coefficients X±(lj) in Eq. ( Bil l. The explicit form of these coefficients are, with Yjo 
\J (21 + l)/47r, expressed as 



X±(w)/Yic 



2V2~K(iQb± 



1 



dr 



r>H\'fs l*^' 



+ ) 



. V2 K qQ rro 



V+ ■ 



dr 



[(t+l)w Jr 



dr 



dr 
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2\Z2nqQ 



in AuT 



dr 



+ 1 



.K 2 fi 



ro 



.dT 



*~S3 I f LL± \ fain AuT 



dr 



in AuT 



— l- 



rfH \ds 
1 d /a±r 



dr 



r 2 drV H $ ±J V ds 1 



AT 



dr , 



(D9) 



for both generic and exceptional modes. In the above expression, V± are the effective potentials of the scalar-type perturbation 
defined in Eq. (|4TT >. and the functions H, Pgi and P53 are as defined in Eqs. ( 148) . The components (t,r) = (T(s),R(s)) 
specifies the trajectory of a charged particle obtained by integrating Eqs. (0|i and <(5j with the initial condition that the particle is 
at rest at the equilibrium position r = r in the limit t — > —00. 

Notice that the integrands in the above expression do not possess any divergences at r = 7'o- By contrast, in the limit r = r + 
the integrands of the last two terms in Eq. (|D9t diverge. To see this, recall that / vanishes at r = r + and dT/ds and d$™/dr 
behave like oc f" 1 there, while the other functions including dR/ds are regular. To ameliorate this singular behavior, we perform 
the integration by parts in the following manner. First, we rewrite e tu)T in the second to last term by using the following identity 



iuiT 



1 dR f de 



itu dT \ dr 



(D10) 



Then, using the relation f (dT/ds) 2 — f 1 (dR/ds) 2 = 1, which is simply the normalization condition of the four velocity, the 
contribution from the last two terms can be rewritten as 



( the last two terms in Eq. 



) 



2k?h 
iuj 

2k 2 /j. 



a ± dT fain 
rH ds ± 



f 



2k> 



+ 

2,, pro 



1- — — ( a±r $ 

r 2 ds dr V H 



dT 
ds 

inAwT 



a±r 

IT 



: dr 



ds 



d (a±r 
~dr \ ~H~ 




$5 



dr 



(Dll) 



where we performed integration by parts in the last equality and dropped the surface terms. The surface terms take the form 
of an infinitely oscillating function such as e lu)T multiplied by a regular factor near the boundaries. Such terms can be safely 
neglected by introducing an infinitesimal damping factor as usual. In the final expression the integrand is free from divergence 
on both boundaries at r = r + and ro. 
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